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Abstract — This paper considers optimal input design wlien 
tlie intended use of the identified model is to construct a 
feed-forward controller based on measurable disturbances. The 
objective is to find a minimum power excitation signal to be used 
in system identification experiment, such that the corresponding 
model-based feed-forward controller guarantees, with a given 
probability, that the variance of the output signal is within given 
specifications. To start with, some low order model problems are 
analytically solved and fundamental properties of the optimal 
input signal solution are presented. The optimal input signal 
contains feed-forward control and depends of the noise model 
and transfer function of the system in a specific way. Next, we 
show how to apply the partial correlation approach to closed 
loop optimal experiment design to the general feed-forward 
problem. A framework for optimal input signal design for feed- 
forward control is presented and numerically evaluated on a 
temperature control problem. 

I. Introduction 

SYSTEM identification is about the estimation and valida- 
tion of mathematical models of dynamical systems from 
experimental data. It is well recognized that the input signal 
used to excite the system during the experiment significantly 
affects the accuracy of the identified model. A correctly 
chosen excitation signal can immensely improve the quality 
of the resulting model while a poorly chosen signal can result 
in a useless model. This motivates input or experiment design 
in system identification. 

The quality of a model depends on the intended use of the 
model. A good experiment should highlight the important 
properties for the intended application. This is the main idea 
of identification for control [1], [2], least costly identification 
[3] and applications oriented input design [4]. 

Much work has been focused on input design when the 
model is to be used in a control application, see for example 
[5] or [6] for application to MPC. However, many industrial 
control systems do not only utilize feedback but also feed- 
forward control. The idea is to measure a disturbance, predict 
its impact on the plant and then compensate for it with the 
input signal. Compensating for measurable disturbances by 
feed-forward control can improve the performance consid- 
erably compared to when only feedback control is used. 
Measurable disturbances could for example be the outside 
temperature when controlling the indoor temperature in a 
house. 
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The aim of this paper is to study optimal input design 
when the intended use of the model is for feed-forward 
control. We will use the application oriented input design 
framework presented in [4]. To illustrate how system prop- 
erties affect the optimal input signal we start by analyzing a 
problem which can be solved analytically. We then present a 
framework for input design for more general systems and 
show how to formulate them as Semi Definite Programs 
(SDP) that can be solved efficiently using numerical methods. 

The outline of the paper is as follows. In Section we 
define the input design problem and give some preliminary 
results from application oriented input design. Section |III] 
analyzes first order FIR subsystems and make some observa- 
tions. A more general framework for input design for feed- 
forward control is given in Section |IV] The framework is 
then applied to a simulation example in Section [V] while 
Section |VT] concludes the paper 

II. Problem Formulation 
Consider the system in Figure [T] The dynamics of the 
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Fig. 1: The considered system. The identification setup 
without feed-forward filter {Ff) and application setup with 
feed-forward filter 

system can be written as 

GQ{q)ut + Hoiq)vt + et, 



Vt 



where q~ is the delay operator, i.e., ut = wt-i- Here 
Go{q) and Ho{q) are the plant transfer functions. The input 
to the system is ut while yt is the measured output disturbed 
by a zero mean white noise signal et with variance Ag. The 
process vt corresponds to a measurable disturbance to the 
system. We assume that vt is a stationary stochastic process 
with known spectral properties, i.e., vt can be written as 
Vt — M {q)st where st is a zero mean Gaussian process with 
variance As and the stable and inversely stable filter M (q) is 
known. The spectrum of vt equals $t,(w) = \M{e^'^)\'^Xs. 

The objective is to design a feed-forward controller Ff in 
Figure [T] to suppress the effect of the measurable disturbance 
in yt. However, the dynamics of the subsystems H and G 
are assumed unknown and need to be identified. Hence we 



want to identify the system dynamics and use the identified 
model to design a feed-forward controller to counteract the 
influence of the measured disturbance vt on the output yt- 
In particular, we will study to design the input ut during the 
experiment such that we with high probability can guarantee 
that the designed feed-forward controller satisfy an accuracy 
constraint. 

We parameterize the submodels as 

G[q,0), Hiq,9), 

where 6 € E," is the model parameter vector and we assume 
that there is a vector do such that 

G{q,9o)^Goiq), H{q,0o) = Ho{q), 

that is, the true system can be described by the model. 

We denote the model parameter vector estimated using a 
Prediction Error Method (PEM), [7], from N data points of 
the inputs and outputs, {ut,vt,yt, t = 1, . . . , N} as 9^. Un- 
der some mild assumptions the asymptotic (in N) covariance 
matrix of the estimated parameters 

P= lim NE\{9r,^9o){9N-9of], 



can be written as 



P- 



1 



27rAe 

where the joint spectrum 



is the spectrum of xo — [^t St] and 



Men^xoi^)j'oie'ndt^, (1) 



T 



dGjq.e) 

de 



M{q) 



♦ dHjqM) 

de 



(2) 



see [7] for details. 

Remark 1: Since we can measure vt, we can correlate the 
input Ut with st (since we know M{q)) and hence $„s could 
be non-zero. 

Further we know that, asymptotically, the identified pa- 
rameter vector 9 Is! lies within the set 



]^[9-9ofp-^[9-9o]<^ 



with probability a. The constant k can be determined from 
the (n) distribution as a function of the probability a and 
n, the number of estimated parameters. 

A. Application oriented input design 

As mentioned earlier we would like to identify the param- 
eters of the system and from this design a feed-forward con- 
troller The feed-forward controller will be denoted Ff{q, 9). 

Using the framework in [4], we let Vapp{9) denote an ap- 
plication cost that measures the degradation in performance 
due to model errors. Here the application cost will represent 
the output variance from the measured disturbance vt when 
designing the feed-forward controller with the identified 



parameters instead of the true parameters, i.e., all acceptable 
models satisfy 

Vapp{,9) ^E{iHo + GoFfi9))vt)^} 



^ r \Ho-GoFj{9)\^<i>,duj<-, 



27r 

where 7 is some positive constant. The requirement is hence 
that the output variance should be less than 1 The overall 
objective of the input design is thus to design our input signal 
to be used during the experiment such that Vapp{9N) < 1/7- 
We can approximate the application cost by a Taylor series 
expansion 

Vapp{e) ^\[o- Oof - ^0] , 

and hence the set of all parameters, 9 that satisfies Vapp (9) < 
1 /7 can be approximated by the ellipsoidal set 



1 



Sapp - {9 : -[9 - 9q] V^'pp [9 - 9o] < - 



(3) 



B. Minimum Variance Input Design 

To satisfy the application constraint we need to insure that 
the identified parameters lies within the application set (|3]) 
with high probability, i.e., that £31 C £app. This is equivalent 
to 



N 



where A ^ B means that ^ — i? is positive semidefinite. 
If the above inequality holds then 9 G Sapp with at least 
probability a. 

The objective here is to find the minimum variance input 
signal used during the system identification experiment such 
that the identified model will satisfy the requirements from 
the application. More formally this can be stated as the 
following optimization problem 



1 

minimize — / <I>,i(a;)(iw 
27r 



N . 



(4) 



subject to -F-i h 7K;p(^o) 

defines a spectrum. 

Assume that we can express the constraint that defines 
a spectrum as an LMl-constraint. Then since in many cases, 
as we will se in the following, P^^ is an affine function of 
the spectral density "I>^„ and since J^^ (l>„(cj)(iaj is linear 
in $^(,, the optimization problem @ becomes a SDP and can 
be solved efficiently. We will come back to how to write the 
spectrum constraint as an LMI later. 

III. First order fir subsystems 

To get some insight in the problem we start by looking at 
a low order example where we can solve the optimization 
problem analytically. We assume that the two subsystems are 
first order FiR-systems, i.e., 

Go{q) = 1 + boq-\ Ho{q) = 1 + /io9~\ 

and the measurable disturbance is zero mean white noise 
with variance i.e., that M — 1 and A„ = A^. We will use 



system identification to identify the parameters — 
and use this to design a feed-forward controller 
The feed-forward controller will be on the form 

( n — t— I 1 f"" ^ 

Ff{q, 



[b h\ 



Co + ciq 



(5) 



1 + flog 1 ' 

where the constants co, ci and ao are functions of h and 
6. Assuming that we know the true underlying system we 
design the optimal feed-forward controller, on the form (0, 
that minimizes the output variance due to the disturbance vt- 
If G{q) is minimum phase then obvious we should invert G 
in the feed-forward filter In the non-minimum phase case it 
is known from minimum variance control that one should 
mirror the zero and then invert the system, see [8]. The 
optimal parametrization is thus given by 

[1 h h] if \h\ < 1 



[co 



ci ao 



hb-'+b- 
P~ 



h 



if \b\ > 1. 



The objective of the system identification is to find an 
estimate 9^ of 6 such that when they are used to design 
a feed-forward controller Ff{9N) the output variance due 
to the disturbance is less than I/7. All parameters, 9, that 
satisfies this is given by 

Vappie) = E{{{H{9o) + G{9o)Ff{e))vtf} < i. (6) 

A. Minimum phase system 

First we look at the case when Go is minimum phase, 
i.e., \b\ < 1. The application cost ^ can now be calculated, 
using for example residue calculus, to 



Vapp{9) 

where c = ho — h + 
then given by 



(P - 2cdb 



A„ < -, 



1-6 7 
b — bo, d = hob — boh. The Hessian is 



2A,; 



P 

-1 



-1 
1 



/?2 



2hobo + 1 



1-bl 



Using ([T]) the asymptotic covariance matrix for the iden- 
tified parameters can be calculated as 



P- 



1 

2^ 



UV 



Idw 



TuiO) r„,(0) 

TuviO) At, 

where r„(fc) = E[utUt-k\ and ruv{k) = E[utVt-k] are the 
autocovariance of ut and the cross covariance between ut 
and Vt, respectively. 

Since the input power can be written as 

1 r 

— J $„(w)da; = r„(0), 

the optimization problem (|4|i can be reformulated as 
minimize r„ (0) 



subject to -- — 

XpK 



P 

-1 



-1 



1 



r„(0) 

ruv (0) Xv 

Note that condition that and r„„ corresponds to a realiz- 
able experiment (or that $^(, defines a spectrum) is that the 



matrix P^^ is positive semidefinite. However, since V^'pp 
is positive semidefinite this is already taken care of in the 
optimization problem. 

The optimization problem above can be solved analyti- 
cally, see [9] for details. The optimization problem above is 
feasible if 

N 

7>0 

AeK 

and the optimal solution is given by 

^XeXvKP 



r„(0) 



N ' 

7AeAi,K 

N ■ 



(7) 



We are now ready to make the following observations; 

• To guarantee that the identified parameters satisfy the 
application requirement with high probability a, we 
require — 7 > 0. Hence the highest possible achiev- 
able accuracy for a given noise variance, experimental 
length and probability is 



7 < 



N 

XpK 



Or conversely, to achieve a certain accuracy the exper- 
imental length need to satisfy 

TV > 7AeK. 

As expected the required input power during the iden- 
tification increases with higher probability (larger k), 
higher noise variance (larger Ae) and tighter application 
requirements (larger 7), while increasing the experimen- 
tal length, N, reduces the required power 
Not as obvious is that the required input power is pro- 
portional to the disturbance variance. One might think 
that a higher variance of Vt could help us to identify 
the parameter corresponding to the disturbance filter, 
H, and hence require less accuracy in the identification 
of G, consequently requiring less power in ut. While 
it is true that a higher power of the disturbance makes 
the identification of H more accurate, the application 
cost is also proportional to Xy. Hence higher accuracy 
is needed to be able to satisfy the application require- 
ments. 

The required input power is proportional to p. The 
constant p captures the effect of the underlying system. 



P 



hi - 2hQbQ 



1 



1-bl 



(ho - bp) 

1 — h'^ 



1 > 1. 



Hence a lower bound on the required input power is 
given by 

r„(0) > 

If the two systems are equal, i.e., bo — ho, then p = 1. 
Thus the least amount of power in the input is needed 
if the two systems are equal. In system identification 
of structured systems it is recognized that it could be 
hard to identify subsystems which are equal, see [10], 



[11]. Here we see that in this particular scenario this 
is not the case. If ho and bo are very different from 
each other p will be large and the required power is 
large. Furthermore we can se that if bo is close to ±1 
then high power is required. This is expected since we 
invert Go in out feed-forward filter and in this case we 
are close to the stability margin and a high accuracy is 
needed. 

If it is known beforehand that the two subsystems are 
equal this should of course be taken into account. In 



this case the feed-forward filter becomes F, 



f 



1 and 



is independent of the identified parameters. Hence it is 
not necessary to perform any system identification. 
The optimal input signal should always be negatively 
correlated with the disturbance and be 

'fXeXvK 



N 



independent of the underlying system. 

B. Non-minimum phase system 

Now we look at the case when Go is a non-minimum phase 
system, i.e., \b\ > 1. The optimal feed-forward controller on 
the form (|5]l will now be 

{hb^ + b-h)q + hb 

^^('^'^) = — W{b^) — • 

Even if we know the parameters of the true system the 
variance will not be zero in this case. The smallest error 
we can make is given by 

{bl-l)ibo-hor 



The application cost is now given by 

Vapp{9)=E[{{H + G,Ff)vtf}^V,r, 



1 

< -. 

7 



Note that I/7 does no longer correspond to the highest 
acceptable variance but is now a bound on how much higher 
the variance can be using the identified parameters compared 
to if we knew the true system. Redoing the same calculations 
as in the minimum phase case, the optimal input signal is if 

b'o 



7- 



> 0, 



given by 



r-n(O) = P2, 



N 



-fXeXvK hoibf, - 3bl + 4) + 26^ - 360 



N 



b'o 



where 



P2 =- 



hlibl 



16bl - 6bf^ - 10) 



+ 



bUbl ~ 1) 
ho{4bl - ml + Ubo) 



4bt 



3bl 



(8) 



b'oibl-i) 

Many of the observations for the minimum phase setting 
also hold in this case, the differences are 



To achieve a certain accuracy the experimental length 
must now satisfy 



b^ 

N > 7AeK-2- 



b'o + ^ 



bt 



The required input power is proportional to p2 defined 
in (|8]l. We see that, for a fixed ho, if |6o| ^ 00 then 



P2 







r„(0) ^0. 



Thus a larger bo requires less input power to achieve a 
certain accuracy. 
• If 60 = ^0 then 

-yXeXvK fep - 6§ + 1 



r„(0) 



N bt ' 

^XeXyK bt-bl + l 



N 



bt 



• Again if 1 then r.u(0) 00. This is at the 

stability margin of our feed-forward controller and it 
is required that we identify G correctly as minimum 
phase or as non-minimum phase. Thus high accuracy 
and consequently high input power is needed. 

C. Signal generation 

We will now show how we can realize a input signal from 
the optimal correlations. Here we will only show how to do 
this for the minimum-phase case as the calculations for the 
non-minimum phase case is analogous. From the optimal 
correlations dTjl we can generate an input signal by using 

ut = Kvt + n, 

where rt is a white zero mean Gaussian process independent 
of vt with variance . The constant K can in this case be 
seen as a feed-forward filter that should be used during the 
identification. Since 

ru{0) = K^X, + Xr, 
ruviO) = KXv, 



we obtain 



K = 



.(0) 



Xr 



Xy 

ru{0) 



X 



M N 



iV2 



>o, 



where the last inequality is due to that N > jX^k to have a 
feasible optimization problem. 

D. White noise input 

Let us compare the results with what we get if we use a 
white input signal during the experiment, uncorrelated with 
the measured disturbance Vf, i.e., r„„(0) = 0. Now the 
optimization problem becomes 



mmimize 



.(0) 



subject to 

XeK 







-1 1 



(9) 



with optimal input 



ru[0) = h 



N 



Comparing to the case with correlated input the required 
power is 

N{N - 7AeK) 

larger, see (|7}. The difference decays as so for large 

N one could instead use an input signal uncorrelated with 

Vt- 

IV. The general case 

We will now study more general cases, when the sub- 
systems are not necessary FiR-filters. Mainly we will focus 
the discussion on how to guarantee that the spectrum we 
optimize over, actually defines a spectrum, and how 
to formulate this as an LMI. Two common ways to do this 
are to use a finite dimensional parametrization or a partial 
correlation parametrization of the spectrum. 

In the finite dimensional parametrization the spectrum is 
written as an infinite series and the optimization parameters 
are given by the truncated vector of coefficients in the series. 
Using the Kalman-Yakubovich-Popov lemma the condition 
that the parameters represent a spectrum can be written as 
an LMI, see [12] for details. 

Often the optimization problem ^ can be written in terms 
of only a finite number of parameters. The idea of the partial 
correlation approach is to find conditions on these parameters 
that guarantee the existence of an infinite extension such that 
the complete sequence defines a spectrum. This can also be 
expressed as an LMI see [12]. 

The two above approaches are however not directly ap- 
plicable to our problem. The problem is that we do not 
have control over the second input Vf. For example, the 
partial correlation approach only guarantees that there exist 
an extension, not that the extension exactly corresponds to 
the given properties of vt- 

The two approaches have however been extended to closed 
loop optimal input design, see for example [13] for the 
finite dimensional parametrization and [14] for the partial 
correlation approach. We will show that our problem is just a 
special case of this. In closed loop input design the objective 
is to find the controller K and the spectrum of the reference 
signal rt, with ut — —Kyt + rt, to achieve some properties 
on the identified models. See Figure |2a] This is the same 
as designing the spectrum of u, and the cross spectrum 
between Ut and St, ^us- This is exactly what we want to do 
in the feed-forward case, we want to find conditions such 
that and corresponds to realizable experiment. We 
see that the feed-forward case in Fi gure l2b] corre sponds to the 
closed loop case when Gq — and Hq — ~M. Hence we 
can use the existing results from closed loop optimal input 
design to guarantee that the spectra in the optimization can 
be realized by a reference signal and a controller K with 
Ut — —Ktjt + rt — Kvt + rt. Here we will show how the 




(a) Feedback 



(b) Feed-forward 



Fig. 2: The relation between the input ut in the closed loop 
case and the feed-forward case. 



partial correlation approach from closed loop optimal input 
design [14] can be used in this setting. 

First we define the generalized moments of the spectrum 



1 

2"^ 



(10) 



, ^^Mland 

de 

is affine in the 



where d{z) = X)/=o ^i'^^ svlcYi that the coefficients are real 
and obey d„i ^ and has all its roots outside the closed unit 
disc. Note that ruk is real and that m_/j ~ mj 

Many common cost functions in input design, for example 
minimum input variance, can be written as linear functions 
of the generalized moments, see [14]. If the model structures 
of G{q, 9) and H{q, 9) are rational then by choosing d{z) 
as the least common denominator of ^'^^^f' 
M*{z) evaluated at 6* = 6*0 we see that P' 
generalized moments too, toi, . . . , to„ c.f., (HJ, see [12] for 
details. 

The input design problem (|4]l can hence be written as a 
SDP in a finite number of generalized moments (fTOl l. Now 
we want to classify the set of all sequences mo,...,TO„ 
that corresponds to a realizable experiment design, i.e., that 
corresponds to a K and $r that are realizable. In [14] a 
semi-definite description of this set is given. We will shortly 
summarize the result here. 

The sequence (mo, mi, . . . , m„) of 2 x 2 matrices satisfy- 
ing TO_fe 



for k — 0, . . . ,n defines a valid experiment 



if the following conditions holcfl 

1) rnk,22 = ^ for k - -n, 

2) YaLo dimk-L2i for A: = 1, . . . , n. 

3) The matrix 



mo 
mi 



TUq 



is positive semi-definite. 
A proof and more details can be found in [14]. 

Remark 2: The above formulation is not only useful for 
feed-forward cases. It could as well be applied, with minor 
modifications, to more general MIMO system identification 
problems where some of the inputs are known but not 



denotes the (p, q) element of the mahix m^.. 



controllable. This is a common case in many industrial 
applications. 

A. Generating the Input Signal 

Solving the optimization problem gives the optimal 
mo, mi, . . . , m„. The problem is then to find the feed- 
forward controller K and the spectrum from theses 
matrices. 

Define the matrix T^{a) = diag(a, 0, a, 0, . . . , a, 0) and 
let Gmax be the largest possible value of a such that 
T„ — T'j[{a) is positive semi-definite. Furthermore let v be 
a non-zero vector such that (T„ — T^(ama£i;))v = with 
V = (p„,q„,p„_i,g„_i,...,po,go)'^- Then one possible 
realization of the reference spectrum and the feed-forward 
controller is 

= |d(e^")pa™a:r, 
K = 

p{z)M{zy 

where p{z) = YIi^oPi^^^ ™d q{z) = I^Lo*^"'- The 
spectrum of r can thus be realizable by filtering white 
noise with variance amax through the filter d*{z). Other 
realizations are possible, see [14] for details and a proof. 

V. Example - House Heating 

In this section we will show how the framework for input 
design for feed-forward control can be applied in a simula- 
tion example. We will consider the problem of controlling the 
temperature in a room using electrical radiators. The control 
signal ut, is the temperature of the radiator and yt is the 
temperature in the room that we would like to control. The 
temperature in the room is also influenced by the outside 
temperature, vt, due to heat transfer in the walls. To simplify 
the problem we consider vt as white noise with variance A„. 
A thermometer is fitted outdoors so we assume that we can 
measure the temperature outside perfectly. The model of the 
system we will use here is a simplified version of the one 
given in [15]. The model is given by 



Vt+i 



'byt + kiUt + k2Vt 



where b, ki and ^2 are some constants that depend on 
the heat transfer coefficients between radiator and air, the 
coefficient through the walls and the sampling time. This 
can be rewritten as 

ki k2 



yt 



l + bq- 



l + bq- 



■ Vt + et 



G H 

where we have added et as zero mean white Gaussian 
measurement noise with variance A^. 

The objective is to identify the parameters 6 — 
(fci,fc2,6) and design a feed-forward filter such that we 
keep the variance in the indoor temperature due to changing 
outdoor temperature less than I/7. 

The feed-forward filter based on the identified parameters 
will be 



1 + bq-^ X'^ + bq 



Since both G and H have the same dynamics we see that 
the accuracy of the the parameter b is not important. The 
application cost is hence 

Vapp{e)=E{{{Ho + GoFf{e))vtf} < i 

which can be calculated analytically using for example 
residue calculus. From this we can then derive V^pp{9o). 

We will now use the framework outlined in Section HV] to 
formulate an optimal experiment for this case. First we look 
at the gradient J'q defined in 



and 



dG(qfi) 



dHjq.e) 

de 





(l+fcq-l)^ 





1 

'{l+bq-^: 



P-1 = 



By defining d{e^'^) = (1 + be^'^Y and m^ as in (fTOl i we can 
write as linear combinations of m^, k = —2, . . . , 2. 
The input signal energy can be written in terms of as 



I r 1 r 



E 

fe=-2 



Skmk,ii, 



-<i>„(w)diLj 



(11) 



where Sk are the coefficients of |c?(e-'")|^, i.e., |d(e-'")p = 
J2k=-2 Ske-'^^ ■ In the same way we can express the output 
variance during the experiment as 



in ^ 
— I ^y{uj)duj ^ ^ l3k [ki fe] mk 



(12) 



where Pk are the coefficients of |1 + 6e^-'"p. 

We can now formulate the optimal input design problem 
where we want to find either a minimum variance input 
signal or a input signal that gives the lowest output vari- 
ance while satisfying the application constraints with high 
probability. The optimization becomes 



mmimize 

mfc,fc=-2,...,2 

subject to 



(Hii) or uni) 



mfc,22 



1 

2^ 



y^^dimk-1,21 
1=0 



mo 
mi 
m2 



mj 
mo 
mi 



|d(eJ")|2 
0, fc = 1, 

>- 0. 



6^ '=", Vfc 



mo 



The above problems are difficult, if not impossible, to solve 
analytically. As the problems are SDPs they can be solved 
efficiently using numerical methods. The following system 
parameters will be used; b — —0.5, A„ = 3, Ae = 1, fci = 1 




Frequency (rad/s) 

Fig. 3: The magnitude of the optimal feed-forward filters to 
be used during the identification for the minimum input vari- 
ance case (i i) and the minimum variance output (i i) . 

For reference Hq is also shown as (i b . 



and ^2 — 0.3. The experimental conditions are = 1000, 
7 = 100 and k = 5.99 corresponding to at least 95% 
probability that the identified model satisfies the application 
requirement. 

The two problems are solved with CVX, a package for 
specifying and solving convex programs [16], giving the 
optimal nik, k = —2, . . . ,2. Using the algorithm outlined in 
Section IIV-AI we generate a reference spectrum and the 
feed-forward controller to be used during the identification, 
K. In both cases the optimal r„ is singular and consequently 
dmax = and <I>r = 0. Therefore it is enough to use 
Ut = Kvt during the identification. The magnitude of the 
optimal feed-forward controllers K for both problems are 
shown in Figure [3] where also Hq is shown for reference. We 
see that when we minimize the input variance the optimal 
K is a lowpass filter, while for the output variance case the 
optimal X is a constant K = —0.3. The case when K = 
—0.3 corresponds to the optimal feed-forward controller 
Ff = — fc2/fci = —0.3 that gives zero output variance in 
y due to the measurable disturbance vt- But using this feed- 
forward filter during the identification gives that yt = et, and 
hence we cannot say anything about the parameters fci, ^2 
and b from this information. However, since we know that 
the subsystems have the same dynamics, this information is 
enough to identify the ratio ^2/^1 which is all information 
we need to design our feed-forward controller! 

We verify the results in a Monte-Carlo simulation. In 
each round system identification data is generated using the 
optimal identification feed-forward filter, K, and the param- 
eters fci, k2 and b are identified. The identified parameters 
are then used to design the feed-forward controller Finally 
we calculate the output variance due to the measurable 
disturbance Vt when the feed-forward filter Ff{9iq) is used 
in the application. In about 98% of the simulations the output 
variance was less than I/7, thus the results seem to be valid. 

VI. Conclusions 

In this paper we have considered optimal input design 
when the identified model will be used in a feed-forward 
control application. First, the first order FiR-filter case was 
considered and some fundamental properties were observed. 
Secondly a framework for optimal input design for feed- 
forward systems was presented and the relation to closed 



loop input design was discussed. The framework was then 
successfully applied to a numerical example. 

Interesting extensions of this work would be to see what 
happens if, on top of the feed-forward, feedback is added 
as this is the most common case in practice. It would also 
be interesting to see if some of the results for the low order 
FiR-case can be extended to more general structures. 

We noted that the framework presented here can be ex- 
tended to more general MIMO system identification problems 
were we only can control a few of the inputs. It would 
be interesting to see if this can be used in identification of 
complex interconnected systems, see [17]. For example if one 
node locally want to identify the dynamics of the network, 
how should it excite the system when it only can measure, 
and not affect, the inputs from its neighbors? 
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